ON GENERALIZED MAX-LINEAR MODELS 
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Abstract. We propose a way how to generate a max-stable process in C[0, 1] 
from a max-stable random vector in R'* by generalizing the max-linear model 
established by Wang and Stoev (2011). It turns out that if the random vector 
follows some finite dimensional distribution of some initial max-stable process, 
the approximating processes converge uniformly to the original process and 
the pointwise mean squared error can be represented in a closed form. The 
obtained results carry over to the case of generalized Pareto processes. The 
introduced method enables the reconstruction of the initial process only from 
a finite set of observation points and, thus, reasonable prediction of max-stable 
processes gets possible. 



1. Introduction and Prelimiaries 

1.1. Introduction. A max-stable process (MSP) ^ = {^t)teK with sample paths 
in C{K) := {/ G : / is continuous} with a compact set if C R has the 
characteristic property that there are for every n g N continuous functions a„ > 
0, hn G C{K) such that 

max (^(*) + 6„)/a„ = ( max {£,t{i) + hn{t)) / a.n{t) ) 

l<i<n yl<i<n J ^^j^ 

where the are independent copies of ^ and "=ci" denotes equality in distribu- 
tion. It is well known (e.g. de Haan and Ferreira (2006)), that MSP are the only 
possible limit processes of linearly standardized maxima of independent and iden- 
tically distributed processes. This is in complete accordance to the well-established 
finite dimensional case of extreme value analysis, see Falk et al. (2010) and again 
de Haan and Ferreira (2006), among others. 

The theory on continuous max-stable processes is essentially based on the early 
works of de Haan (1984), Gine et al. (1990) and de Haan and Lin (2001), followed 
by recent findings with the focus on different aspects within the theory, e.g. Hult 
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and Lindskog (2005, 2006), Stoev and Taqqu (2005), Davis and Mikosch (2008), 
Kabluchko (2009), Wang and Stocv (2010), Aulbach ct al. (2012). 

Moreover, the class of excursion stable generalized Parcto processes, which is 
closely related to the class of max-stable processes, was examined in Buishand 
et al. (2008) and Aulbach et al. (2012), the most recent advances on that issue can 
be found in Ferreira and de Haan (2012). 

There is a very crucial problem of the theory on stochastic processes concerning 
its relevance in practice: as (continuous) processes as a whole cannot be measured 
exactly, the question arises how to construct those processes (with some character- 
istic stochastic behavior such as max-stability) from a finite set of observations. As 
an example one can think of data from a finite set of measuring stations measuring 
the sea level along a coast and one is interested in predicting the sea level between 
those measuring stations. 

This is an issue of conditional sampling and in the case of max-stable processes 
(or fields, if the domain has more than one dimension) there are (partial) answers 
on the arising questions in Wang and Stoev (2011) and Dombry et al. (2012). 

In this paper, we pick up the so-called "max-linear model" introduced in Wang 
and Stoev (2011): for appropriately chosen nonnegativc continous functions 
(7o, ■ • ■ , 3d one obtains an MSP r] ~ {'rit)te[o,i] by setting 

r]s = max s £ [0, 1], 

where X = {Xq, . . . , Xd) is a max-stable rv with independent components. The 
obvious restriction of this model is the required independence of the margins of X 
which results in the fact that rj has always a discrete spectral measure (the latter is 
in our setup descripted below essentially the distribution of the generator process 
Z; see Aulbach et al. (2012) for details). 

In Section 2, we generalize this model by allowing arbitrary dependence struc- 
tures of the margins of the max-stable rv X. This leads immediately to the main 
issue of this paper, the reconstruction of a max-stable process which is observed 
only through a finite set of indices: it is shown in Section 3 that if the random vector 
is some finite dimensional projection of some initial MSP, the processes resulting 
from the construction in Section 2 converge uniformly to the original process as the 
grid of indices gets finer. Moreover, the mean squared error between the predictive 
and the original process is computed at a fixed index, which is useful for practical 
purposes. 

It is also possible to predict generalized Pareto processes with the same tech- 
niques as for MSP which is the content of Section 4. 
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To begin with, the following subsection recalls some basic theory needed in what 
follows and introduces some notation. For the ease of notation we choose K = [0, 1] 
as domain of the processes, being aware of the fact that all results arc valid for an 
arbitrary compact set A' C M as well. The extension of the results to more general 
domains (in particular higher dimensions of the domain) is not immediately obvious 
and subject of current research. 

1.2. Prelimiaries. We call a random vector (rv) X = {XQ,...,Xd) standard 
max-stable, if it is max-stable and each component follows the standard nega- 
tive exponential distribution, i.e., P{Xi < x) = exp(x), x < 0, i ~ 0, . . . , d. It is 
well-known (e.g. de Haan and Resnick (1977), Pickands (1981), Falk et al. (2010, 
Sections 4.2, 4.3)) that X is standard max-stable if and only if (iff) there exists a 
rv Z = (Zq, . . . , Zd) with Zi £ [0, c] almost surely and E{Zi) = 1, i = 0, . . . , d, for 
some number c > 1, such that 



P{X <x)^ cxp(- ||a;||^) := exp -E max (|x,| Zi)]] , a; < e K'' 



Note that ||-||^ defines a norm on M'^'+^. called ZJ-norm, with generator Z. The D 
means dependence: We have independence of the margins of X iff ||-||p equals the 
norm ||a;||j^ = X]f=o I'*"*!' which is generated by (Zq, . . . , Zd) being a random permu- 
tation of the vector (d -I- 1, . . . , 0). We have complete dependence of the margins 
of X iff is the maximum-norm ||a;||oo ~ maxo<i<d which is generated by 
the constant vector (Zq . . . , Zd) — (1, . . . , 1). We refer to Falk et al. (2010, Section 
4.4) for further details of D-norms. 

We call a stochastic process with sample paths in ^^[0,1] := {/ G C[0, 1] : 
/ < 0} a standard max-stable process (SMSP) , if it is a max-stable process with 
standard negative exponential univariate margins. 

Denote by E[Q, 1] the set of those bounded functions / : [0, 1] — !> R that have only 
a finite number of discontinuities, and let -E~[0, 1] be the subset of those functions 
in -EfO, 1], which attain only non-positive values. 

It is known from Gine et al. (1990) and Aulbach et al. (2012, Lemma 2) that 
a stochastic process rj in C[0, 1] is an SMSP iff there exists a stochastic process 
Z = (Zf)4g:[o,i] with sample paths in C'+[0, 1] := {/ G C[0, 1] : / > 0} with Zt < c 
a.s. and E{Zt) = 1, t G [0, 1], for some constant c > 1, such that 



P{ri < /) = cxp(- 11/11^) := exp -E sup {\f{t)\ Z*) , / G E-[0, 1]. 







The condition P(supjg[Q ij Zt < c) = 1 on the generator process Z can be weakened 
to E ( supt£[o,i] ^t] < 00, see de Haan and Ferreira (2006, CoroUary 9.4.5). 
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Note that ||-||^ defines a norm again, this time on the space i?[0,l]. It is also 
called £)-norm with generator process Z ^ and we have 

||/IL<ll/lb<£i^ll/lloo, /ei?[0,i], 

where ||/||^ = suptg[o,i] l/WI = = ^(ll^lloo) is the extremal 

coefficient , cf. Smith (1990). This implies |H|^ = |H|^ iff Sd = 1, cf. Aulbach 
et al. (2012). Moreover, the preceding inequality shows that each D-norm on the 
space -E[0, 1] is equivalent to the sup-norm || ■ || , which is itself a D-norm by putting 
Zt = 1, t G [0,1]. 

We conclude this section by introducing generalized Pareto processe as consid- 
ered in Section 4. For the purpose of this paper, the following definition is suf- 
ficient: we call a stochastic process V in (7" [0,1] a standard generalized Pareto 
process (SGPP), if there exists a D-norm ||-||^ on i?[0, 1] and some c > 0, such 
that P(y < /) = 1 - 11/11^ for all / e i?-[0,l] with ||/||^ < c. For a detailed 
examination of GPP we refer to Ferreira and de Haan (2012). 

2. The generalized max-linear model 

Let X = (Xo,...,Xc;) be a standard max-stable rv with pertaining D-norm 
||-||^^ ^ on K''+i generated by Z = (Zq, . . . , Z^), d e N, i.e. 



P (X < a;) = cxp 1^— ^ j = cxp \^ max |a;i| 

X = [xo, ■ . ■ , Xd) < 0. Choose deterministic functions go, ■ ■ ■ , gj, ^ [0, 1] with the 
property 

(1) ||(5o(t),...,gd(t))|lz5„,...,, = 1, t6[0,l]. 

For instance, in case of independent margins of we have ^ = ll'llii and 

condition (1) becomes 



Y,9^{t) = l, tG[0,l], 



i=0 



i. e. gi{t), i = 0, . . . , d, defines a probability distribution on the set {0, . . . , c?} for 
each t e [0, 1]. In this case, an example is given by the binomial distribution 

(%\l~tf-\ i = Q,...,d, te[0,l]. 



Put now for t £ [0, 1] 

(2) Tjt := max 



»=o,...,d gi{t) 

The model (2) is called generalized max-linear model . It defines an SMSP as the 
next lemma shows: 
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Lemma 2.1. The stochastic process rj = (?7t)te[o.i] ^'^ (2) defines an SMSP with 
generator process Z — (Zt)t<z[o,i] given by 

Zt = max ig^it)Z,), t e [0,1]. 

i—0,...,d 

Proof. At first wc verify that tire process Z is indeed a generator process. It is 
obvious tliat tfie sample patlis of Z are in (7+ [0,1]. Furthermore, we have by 
construction for each t G [0, 1] 



EiZt)^\\igoit),...,gdm\D,,. 



= 1. 



As 



< 



ID 



for an arbitrary D-norm, we have ||(go(i), ■ ■ ■ t 9d{t))\\^ < 1, t 6 



[0, 1], and, thus, Zt < maxi=o,....d Zi, t G [0, 1]. 
In addition, we have for / 6 E~ [0, 1] 

PiV<f)^P{X^<g^it)fit), Z=0,...,d, is [0,1]) 



= P[X,< inf i^O,...,d 

*e[o,i] 



P{X,<- sup i-O,. 

V *6[0,1] 



exp 



sup iga\fit)\),..., sup {gd{t)\fm 
,te[o.i] ie[o,i] , 



Do 



exp -i; max sup {gi{t) \f{t)\) Zi 

\ y=o,...,d \ tg[o,i] 



exp -E sup \f{t)\ max {gi{t)Zi) 



exp sup {\f{t)\Zt 

V Vte[o,i 



which completes the proof. 



□ 



Remark 2.2. Condition (1) ensures that the univariate margins rjt, t G [0,1], of 
the process rj in model (2) follow the standard negative exponential distribution 
P{ilt 1^ x) — exp(x), X < 0. If we drop this condition, we still obtain a max-stable 
process: Take for n G N i. i. d. copies 77^^^ . . . , tj'"' of ij. We have for / G E~ [0, 1] 



P (^n max ■q^^'^ < = P [r] < - 



P[Xi< inf 

te[o,il 



i = 0, 
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= exp ^- 

= exp ^- 

= P(V < f) 
The univariate margins of r) are now given by 

(3) P(77t<a;)=cxp(||(,9o(i),---,.9dW)llDo,. .,. ••^)' 2; < 0, te[0,l]. 

Note that the above calculations also give an alternative proof of Lemma 2.1, except 
that wc do not obtain the generator process of t] with this approach. 

In model (2) we have not made any further assumptions on the D-norm ^, 
that is, on the dependence structure of the random variables Xq, . . . , X4. The spe- 
cial case W-Wug ^ = IMIi characterizes the case where Xq, . . . ,Xd are independent. 
This is the regular max-linear model , cf. Wang and Stoev (2011). 

On the contrary, ^ = II '11 00 provides the case of complete dependence 

Xq = ■■■ = Xd a. s. with the constant generator Zq = ■■■ = Zd = 1. Thus, 
condition (1) becomes maxi=o,...,d ~ 1, t £ [0, 1], and therefore 

Zt ^ max {gi{t)Zi)^ max ^^(t) = 1, ie[0,l]. 

0,...,d 0,...,d 



sup 

vte[04] 



5o|./(t)| 



, sup 

te[o,i] 



9d{t)\f{t)\ 



Do d. 



sup (ffo |/(i)|) , ■ • ■ , sup igdit)\fit)\ 
te[04] te[o,i] 



-Dn 



3. Reconstruction of SMSP 

The preceding approach offers a way to reconstruct an SMSP in an appropriate 
way. Let r) = (?7t)tG[o,i] be an SMSP with generator process Z = {Zt)t^[t)^i] and 
D -norm Choose a grid — : sq < si < • • ■ < Sd-i <■ Sd '■ — 1 of points within 

[0,1]. Then (t/sq , . . . , t^s^ ) is a standard max-stable rv in M'*+^ with pertaining 
D-norm IMI^,,, ^ generated by {Zs„, . . . , ZsJ. 

The aim of this section is to define some SMSP i) = (f/t)te[o.i] for which fis^ ~ 
Vsi, i = 0,...,(i, holds, i.e. 17 interpolates the finite dimensional projections 
{rjsa , . . . , rjs^ ) of the original SMSP t] in an appropriate way. This will be done 
by means of a special case of the generalized max-linear model and we show that 
this way of "reconstructing" the original MSP rj is reasonable as the pointwise 
mean squared error MSE (jjf'^^^ := E (^(j]t — Vt'^^^ ^ diminishes for all t G [0, 1] as 
d increases. Indeed, there is additionally uniform convergence of the "predictive" 
processes and the corresponding generator processes to the original ones. 
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3.1. Uniform convergence of the discretized versions. As we have shown 
in Lemma 2.1, the stochastic process t) = (?70te[o,i] in (2) defines an SMSP with 
generator process Z = (Zf)fg[o,i]> given by 

Zt - max {g,{t)Z,^), te [0,1]. 

i—0,...,d 

Denote by ||-||^. ^ . the _D-norm pertaining to the bivariate rv {r]s-_-^,r]s-), i = 
1, . . . , d. Define tlie functions go, ■ ■ ■ , gd by 

Sl-t 



50 (i) := 



(Sl ~t,t)\\jj^^^' 



0, 



t e [o,si], 

else. 



t — S,_l 



Si+l - t 



i+l -t,t- St)\\jj^ .^^ ' 



0, 



9dit) := 



{Sd-t,t- Sd-l)\\o^_^^ 



0, 



t e [si_i. Si], 

ie[s,,s,+i], i = l,...,d-l, 
else, 

t G [sd-i, Sd], 

else. 



It is obvious that go,. . . ,gd G ^+[0, 1]. Moreover, we have for t S [si-i,Si], i — 
l,...,d, 

Ugoit), . . • ,.9dW)llA,,...,. = = 1- 

Hence, the functions go, ■ ■ ■ ,gd ai'c suitable for the generalized max-lincar model 
(2). In addition, they have the following property: 

Lemma 3.1. The functions go, . . . ,gd defined above satisfy 

llfilloo = S^sO = 1, i = 0,...,d. 

Proof. From the fact that a Z?-norm is monotone and standardized we obtain for 
i = 1, . . . , d — 1 and t G [si-i, Si) 

t-s,^i 1 , 1 



\\{si ~t,t- Si- 



,1 



< 



-D,-i, 



ll(o,i)lli.._,. 



1, 



and for i G [si, Si+i) 



g^it) 



Si+1 - t 



[St+1 



t,t ~ s. 



< 



l(l,< 



Di 



Analogously, we have go < 1 and ffd < 1. The assertion now follows since gi{si) = 1, 
i = 0,...,d. □ 
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The SMSP t) = ('7t)te[o.i] that is generated by the generahzed max-hnear model 
with these particular functions go, . . . ,gd is given by 

(4) 



rjt = max 



= II (sj ~t,t - Si-i)\\jj^_^ ^ max 1 ""2^ -rz^ — ) > ^ G [sj-i, Si], i = 1, . . . , d. 

Note that r\s^ < almost sure, i ~ 0, . . . , d. This implies that the maximum taken 
over d + 1 points in (2) goes down to a maximum taken over only two points in (4) 
since all except two of the gi vanish in t G [si-i, Si], i ~ I, . . . ,d. 

The above process obviously interpolates the rv {rjg^ , ■ • ■ , ?7srf ) ■ In summary, we 
have proven the following result. 

Corollary 3.2. Let r] = (?7t)tG[o,i] be an SMSP with generator Z = (^t)te[o.i] ,• 
and let :^ sq < si <,...,< s^-i < Sd '■= I be a grid of points in the interval 
[0, 1]. The process i) ~ (?)t)tG[o,i] defined in (4) is an SMSP with generator process 
Z = (Zt)tg[o,i], where 

(5) Zt =iTL&x{^gi^i{t)Zs,_i,gi{t)Zs^, t € [si^i, Si], i = 1, . . . ,d. 

The processes i) and Z interpolate the rv (r]sg, . . . ,r]s^) and {Zs^, ■ ■ ■ ,Zs^), respec- 
tively. 

We call fj the discretized version of ij and Z the discretized version of Z, both 
with grid {so, ■ • • ,Sd}- We show now that the preceding approach is suitable to 
approximate the underlying SMSP solely by multivariate observations; that is, the 
discretized version of the underlying SMSP converges to this very process in a strong 
sense. We need the following two lemmata which provide some technical insight in 
the strucure of the chosen max-linear model. 

Lemma 3.3. The SMSP defined in (4) fulfills fori = l,...,d 

sup f)( = max (?7s,_i,77si) , 
te[si_i,si] 



inf ^fjt ^ -\\{r]s,_,,r]s,)\\^ 
This minimum is attained for t = (si_i77s._j + Sirjs^) / {r]s^_^ +Vsi)- 

Proof. We know from Lemma 3.1 that gi^i{t), gi{t) < 1 for an arbitrary i = 1, . . . , d 
and t € [si_i,s,;]. Hence, 

^t = max -— , — -— < max(77,,,_i,?7^J 

\gi-i{t) gi{t)J 
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for i = 1, . . . , c? and t G [si-i, Si], which yields the first part of the assertion. Recall 
that j]si < with probability one, i = . . . ,d. 
Moreover, we have for t G {si-i, Si) 



^ < 



t > 



where equality in one of these expressions occurs iff it does in the other two. In 
this case of equality we have 

'7s 



t, t S-i — 



\i'ns,-i,Vs,)\ 



t - Si^l 

On the other hand, the nionotonicity of a Z?-norm implies for every t £ {si-i,Si) 
with t > (si_i77s^_i + s^r]s^)/ii]si_i + Vs^) 



fjt = \\iSi~t,t- S.,-i)\\jj^_^^^ ■ 
Si - t 



t - Si 



> 



t - Sj_l 



,1 



Di- 



^,1 



Vs 



Di- 



Recall again that ry^. < almost sure. The case t < {si-iris-_-^ +Si?7si)/(^si_i +?7si) 
works analogously. □ 

As an immediate consequence of the preceding result we obtain that for x < 

i) <x max (Tys^, . . . , rjs^) < x 

and 

rj > X <^=^ max ||(77si_i,'7s,)||n < -a::- 
The next lemma is on the structure of the underlying generator processes. 

Lemma 3.4. The generator process defined in (5) fulfills for i ~ I, . . . ,d 

sup Zt = max {Zs^_-^ ,Zs.) . 
te[si_i,si] 



In particular, the extremal coefficient E 



of the SMSP i) coincides with 



the extremal coefficient E(iiiajii^o,....d Zsi) of the rv (r/s^, . . . , 77^^). Moreover, for 
i = l,...,d, 

inf /("^^^^''-^'^^^''^"^'--)"' z/ Z,, >0, 

* \o else. 

In the first case, the minimum is attained for t = (si_iZs.+SiZs._j)/ {Zg^_^ + Zs^- 
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Proof. We know from Lemma 3.1 that gi-i{t), gi{t) < 1 for an arbitrary i = 1, . . . , d 
and t € [si_i,Si]. Hence, 



{gi^i{t)Zsi_^,gi{t)Zs^ < max(Zs^_i,Zs,) 



for i ~ 1, . . . ,d and t G Si], which yields the first part of the assertion. 

Moreover, we have for t G (si_i, s,) in case of Zg._-^, Zg. > 

Si-lZg- + SiZg. -^^ 



{si - <it~ s,_i)Zs 



Sj — t ^ Zg 



t — Si- 



t > 



Zg 



where equahty in one of these expressions occurs iff it does in the other two. In 
this case of equahty we have 



On the other hand, the monotonicity of a Z?-norm implies for every t G (si_i,Si) 
with t > (.s,_iZ,, + s,Zg^_,)/{Zg^_, + Zg^) 



Si - < 



> 



t - Sj-l 



1 



Di- 



\il/Zg^_„l/Zg,)\\^^_^^ 



The case t < {si^iZg. + SiZs._^) / [Zg^ ^ + Zg.) is shown analogously. □ 

So far we have only considered a fixed discretized version of an SMSP. The 
next step is to examine a sequence of discretized versions with certain grids whose 
fineness converges to zero. It turns out that such a sequence converges to the initial 
SMSP in the function space C[0, 1] equipped with the sup-norm. Thus, our method 
is suitable to reconstruct the initial process. 

Let 

g,:={s['\s['K..^sf}^ 0=:.r<4')<...<4'):=l, d G N, 
be a sequence of grids in [0, 1] with fineness 

>d^oo 0. 



, (d) (rf) 
max \s) ' — s)_^ 
= l,...,£i ^ ' ' ^ 
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Let fj^''^ = (?7t''^)tg[o.i] be the discretized version of an SMSP rj — (?7t)te[o,i] with 
grid Qd- Denote by Z^'') = (^t''^)tg[o,i] and Z = (^t)tG[o,i] ^^e generator processes 
pertaining to f)'''-' and J7, respectively. 

Theorem 3.5. The processes i)^'^^ and Z^'^\ d G N, converge uniformly to ij and 
Z, respectively, i. e. ||i7^''^ ~ ''?|| ~^d->oo and Z^''-^ — Z 0. 

Proof. Denote hy [t]d, d £ N, the left neighbor of t £ [0, 1] among Qd, and by {t)d, 
d G N, the right neighbor of t G [0, 1] among Qd- Choose a sequence s^''^ G [0, 1], d G 
N, with s('') -^dGN s G [0, 1]. Then obviously [s'^'^^]d ^d^oo s and {s'^'^^)d -^d^oo s. 
Hence we obtain by Lemma 3.3, and the continuity of the process rj 

se[[s<''']d,(s('''>d] 

as well as 

which yields the first part of the assertion. 

Now we show that Z*'*' -^d^oo Z in (C[0, 1], IMI^o)- K ^ 0, the continuity 
of Z implies Z^g(d}]^ 7^ 7^ Z^^fd)^^ for sufficiently large values of d. Repeating the 
above arguments, the assertion now follows by Lemma 3.4. If Zg = 0, the continuity 
of Z implies 

^sfi) < 2max(Z[g(d)]^,Z|;^(d)^^) ^d^oo = 0, 

which completes the proof. Check that ||((s'-'^'')d ^ ^ [s'-'''']fi)||j) > 1/2 since 
every Z?-norm is monotone and standardized. □ 

The preceding theorem is the main reason why we consider the discretized version 
r) of an SMSP t] a reasonable predictor of this process. The predictions fjt of the 
points rjt, t G [0,1], only depend on the multivariate observations {rjsg, . . . ,r]s^). 
More precisely, the only thing we need to know to make these predictions is the 
set of the adjacent bivariate marginal distributions of (rjsg, . ■ . ,r]s^), that is, the 
bivariate D- norms |H|£i._j ., i = l,...,d. Since it is assumed that the values t]s^, 
i = 0, . . . ,d are known, the prediction fit usefully coincides with the actual value rjt 
for t G {so, ■ • ■,Sd}. 

3.2. The Mean Squared Error of the Discretized Version. Considering % a 
predictor of Tyt we can ask for further properties such as the mean squared error, 
which is our next aim. For that purpose, we formulate the following lemma. It 
applies to bivariate standard max-stable rv in general. 
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Lemma 3.6. Let [X, Y) he a bivariate standard max-stahle rv, i. e. there exists 
some D-norm such that P{X < x,Y <y) ~ cxp (— || (a;, y)!!^)), a;, y < 0. Then 



E{XY)^ / 2 



dt. 



In particular, the covariance and the correlation coefficient g of X and Y are given 
by 

Cov(X, Y) = r dt-1^ g{X, Y). 

Proof. Recall that the expected value of a negative exponentially distributed ran- 
dom variable ^ with parameter A > is given by 



I a;Aexp(Ax) dx = — 1/A. 

J —oo 



Hence, elementary calculations show 

/o /-O 
/ P{X <x,Y <y) dx dy 
-OQ J —OQ 

M 

I cxp(- \\{:x,y)\\jj) dx dy 

^0 

/ exp(a;||(l,y/a;)||^) dx dy 
x exp (a; 11(1, dx du 



1 





a; II (1, u) lljj cxp (x II (1, u) 11^) dx du 



du. 



II(i,^)IId 



The remaining assertions follow from the fact that E(X) = E{Y) = —1 and 
Var(X) = Var(r) ==1. □ 

Example 3.7. In accordance to the characterization of the independence and com- 
plete dependence case in terms of _D-norms, we obtain in case of ||-||^ = \\-\\i 



Co.iX,Y)=f^^du-l = 



and in case of \\-\\]j 



Cov{X,Y)= 2^^-1 = 1. 

Jq (max(it, 1)) 

In particular, we have Cov{X,Y) ~ g{X,Y) e [0,1] for every bivariate standard 
max-stable rv {X, Y) since the maximum norm is the least D-norm and the sum 
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norm is the largest D-norm. In addition to this, we obtain for 



ID — 11112 



f°° 1 r 1°° 

Cov{X, Y)^ — du-l= arctan(it) - 1 = 7r/2 - 1. 

Jq [u +1) I 







In order to calculate the mean squared error of the predictor fit, we have to 
determine the mixed moment E{'r]tf}t)- For this purpose, we want to apply the 
previous lemma which is why we need to ensure that the vector (774, f]t) is standard 
max-stablc itself. This is the content of the following Lemma. 

Lemma 3.8. Let t] ~ ('/t)te[o,i] ^6 SMSP and denote by rj — (jit)te[o,i] 
discretized version with grid {sq, . . . , .s^}. Then the bivariate rv (r]t, fit) is standard 
max-stable for every t E [0, 1] with the pertaining D-norm 

\\{x^y)\\Dt-= {x,9i-i{t)y,9i{t)y] , i e [si_i,Si], i = 

^ ^ DtA-l.i 

where IMI^jj . ^ . is the D-norm pertaining to (jitT^si^nVsi)- 

Proof. We have for every t G [si-i, Si], x,y < and i = 1, . . . , d 

P{Vt < x,T]t <y)^ P{vt < x,r]si_, < gi-i{t)y,i]si < gt{t)y) 

= exp ( - e( max ( \x\ Zt,gi-i{t) \y\ Zs^_-,,gi{t) \y\ Z^, 



= exp - E(^ma,x \x\ Zt, \y\ max {gi^i{t)Zs^_i,gi{t)Zs. 

The vector 

{^Zt , max 1 (t) _ J , Qi {t) Zsi)^ 
defines a generator for every t € [si_i, s;], i = 1, . . . , d as for all such t 

;(^max(5i_i(i)Zs._i,5i(t)Zs.)j = IKSi-i = 1- 



E 



□ 



Let us recall the sequence of processes we have discussed in Theorem 3.5. Sup- 
pose ri is an SMSP and choose a sequence of grids Qd of the interval [0, 1] with 
fineness Kd — J-d^-oo 0. Denote by d G N, the sequence of dicretized versions 
of 77 with grid Qd- Denote further by |H|j-,(d) the Z?-norm pertaining to {r]t,fjl'^^), 
te [0,1], d£N. 

Theorem 3.9. Let r] and fi^'^\ d G N, fee as above. The mean squared error of fj^^ 
is given by 

MSE (r)^^) := E ( (,t - ) = 2 | 2 - / ; du \ ^d- 
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Proof. The second moment of a standard negative exponentially distributed random 
variable is two and, therefore, we obtain by Lemma 3.6 and Lemma 3.8 



= 4 - 2 / 2 — 

"'O ||(l,lt)||^(d) 

Next, we show IHI^(d) ->d^oo ll-lloo for all t G [0, 1]. Denote by Z and Z^, d e N, 
the generator processes of r) and r)'''-', c? G N. Define 

P / y\ ^ ,7 suPtg[o.i] 
m :^ h \ sup Zt \ < oo and Z := — ^ . 

V*6[o,i] / m 

Then E{Z) ~ 1 and, thus, {Zt, Z) defines a generator for all t G [0, 1]. Denote by 
||-||^ the £)-iiorm pertaining to this generator. Lemma 3.1 and Corollary 3.2 imply 
^(t^) < z for all d G N. Therefore, we have for arbitrary x, G M, c? G N and 
t G [0, 1] 

max {^x\ Zt, \y\ Z^^^ < max (^\x\ Zt, \my\ Zt^ , 

where 

E (max ( |a;| Zt, \my\ Zt]) = \\{x,my)\\Q < oo. 



Hence, the dominated convergence theorem, together with the fact that Z^'*-' — ^^-s-k 
Zt for all t G [0,1] by Theorem 3.5, implies 

\\{x, = E (max Zt, \y\ ^^^^o E (max(|a;| Z*, \y\ Zt)) = ||(a;,2/)||^ 

for all x,y gR. 

In Example 3.7, we have already calculated IKlj'")!!^^ = 2. Since \\-\\^ 
is the least £)-norm, we have for all d G N and t G [0, 1] 

1 1 

< 



||(l,u)||^(., - \\il,u)fj 

and therefore by the dominated convergence theorem again 

/■°° 1 f°° 1 

^0 ||(l,^)f^(., io |l(l,u)||L 

which completes the proof. □ 
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4. Reconstruction of SGPP 

The preceding technique of discretizing and reconstructing a given SMSP can 
also be applied to the case of SGPP by simply replacing the standard max-stable 
rv in the model (2) by a standard generalized Pareto distributed rv. Again, the 
generalized max-linear model results in an SGPP. Once this statement is proven, 
most of the results of the previous sections carry over in a very straightforward way. 

Recall that a stochastic process V in ^"[0, 1] is an SGPP, if there exists a D- 
norm ||-||^ on i?[0, 1] and some cq > 0, such that P{V <./) = ! — WfWo ^o'' 
feE-[0,l] with ||/|L<co. 

4.1. Uniform Convergence of the Discretized Versions. Let {Yq, Yi, . . . , Yd) 
follow a standard generalized Pareto distribution (GPD), i.e. there is a ZJ-norm 

IMI A, , on K'^+i generated by (Zo, . . . , Z^), and a vector = (y^"', . . . , yf^) < 

0, such that P{Yo < yo, ■ • • , >d < 2/d) = 1 - llylli,^ ^ for aU y = {yo, ...,yd) with 
y(o) < y < 0. Note that this implies that each univariate marginal distribution 
of {Yq, . . . ,Yd) coincides in the upper tail with the uniform distribution on [—1,0]. 
For a detailed examination of GPD rv, see e.g. Falk ct al. (2010). 

We apply the generalized max-linear model to the GPD rv and obtain a stochas- 
tic process V = {Vt)te[o,i], 

Yt 

(6) Vt := max — 

i=0....,d gi[t) 

where go, - ■ ■ ,gd G C'^P, 1] are functions satisfying condition (1). Again, this pro- 
cess defines an SGPP as the next lemma shows. 

Lemma 4.1. The stochastic process V ~ {Vt)te[o,i] *^ (6) defines an SGPP with 
generator process Z = [Zt)te[o,i]} 

Zt = max ig.,it)Z,), t £ [0,1]. 

i—Q,...,d 

Proof. We have already shown in Lemma 2.1 that Z defines a generator process. 
Put Co := - max,=o,...,d (j/f V llffjiloo) • Then we have for ah / e .^"[0, 1] 

II./IL < ^0 ^ ^ inf^j f{t) > ^. max^ [yf/ \\g,\Q 
and therefore for i = 0, . . . , d 

inf ^ (.9,(i)/(i)) > inf , { g^{t) . max _ f \ 



te[o,i] ■ te[o,i] \ i=^'---'d\\\gf 

( infte[o,i]5»(i) (o)^ 

= max — - • y) 

j=o....,d \ suptg[o_i]5j(t) J 
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^ inffg[o,i] 9i{t) ^ (0) 
~ suptg[o,i]5iW ' 

^ (0) 

Hence, we have for all / close enough to zero 

P{V <f)^P{Y,< g.{t)f{t), * - 0, . . . , d, t e [0, 1 



P[Y,< inf (5,(t)/(t)), * = 0,.. 
te[o,i] 



1 - 



sup (<?o 1/(01),- sup {gd{t)\f{t)\) 
vte[o,il te[o,i] 



l-E\ max sup |/(t)|) 
\i=o,...,d \ tg[oa] 



1-S sup 1/(01 max (5,(t)Z,) 
\te[o,i] V 

sup (1/(01 ZO) , 



which completes the proof. 



□ 



Just like ill the SMSP case, we can use this model to discretize and reconstruct a 
given SGPP. Let V = {Vt)t^[o,i\ be an SGPP with generator process Z = (Zt)tg[o^i] 
and D-norm Choose a grid =: so < si < • ■ • < Sd-i < Sd := 1 of points 

within [0, 1]. Then {Vsg, . . . , Vg^) is a standard GPD rv in M''+^ with pertaining 



D-norm 11 -I 



generated by (Z^, 



Now choose deterministic functions go, ■ ■ ■ ,gd G C^[0, 1] with the property (1) 
and put for t £ [0, 1] 

v.. 



(7) Vt := m.... . 

i=0,...M gi(t) 

As we have shown in Lemma 4.1, the stochastic process V = (V4)te[o,i] in (7) 
defines an SGPP with generator process Z ~ {Zt)t£io.i]: 

Zt = max (5,(0^s.), ^ e [0,1]. 

i—0,...,d 

By choosing the exact same functions go, . . . ,gd as in the special case of the gener- 
alized max-linear model in Section 3, we obtain the process 
(7') 



Vf = max 



.9i-i(i)' 5i(0 
II (sj ~t,t - Si_i)||^^_^ ^ max 



Vs. 



Vs. 



t t-s 



i-l 



t G [sj_i,Si], i = 1,.. 



ON GENERALIZED MAX-LINEAR MODELS 



17 



In order to show that this process defines an SGPP as well, we only have to verify 
that the functions go, . . . ,gd realize in C''''[0, 1] and satisfy condition (1), which we 
have already done in Section 3. Thus, the following result is proven. 

Corollary 4.2. Let V = (?7t)tG[o,i] SGPP with generator Z ~ {Zt)te[o.i], o-i^d 

:— So < si <,..., < Sd-i < Sd '.^ I be a grid in the interval [0, 1]. The process 

V = {Vt)te[o,i] defined in (7') is an SGPP with generator process Z ~ {Zt)te[o.i]; 
where 

Zt ^ ma,x{gi_i{t)Zsi_i,gi{t)Zsi) , ie[sj_i,Si], i d. 

The processes V and Z interpolate the rv (T4o, • ■ • , Vs^) and {Zsg, . • • , Zs^), respec- 
tively. 

In complete accordance to the SMSP case we call V the discretized version of 

V with grid {sq, . . . ,Sd}. 

Remark 4.3. Let the original SGPP V satisfy P{V </) = !- ||/||^ for all 
/ G E~[0, 1] with ll/lloo < Co, with some D-norm ||-||^ and some cq > 0. It is clear 
that this implies 

PiVso<yo,-.-,Vs,<yd)^l-\\iyo,...,yd)\\Do,„„, 

for all y := {yo, ...,yd) with -cq < < 0, z = 0, . . . , d. Now denote by the 
ZJ-norm pertaining to the discretized version V of V with grid {sq, . . . , Sd}- Just 
like in the proof of Lemma 4.1 we obtain P{V < /) = 1 — for all / G E^[0, 1] 

with 11/11^ < Co, where 

Co = - max - — — = Co, 
*=o.-.All.9dloo/ 

since in this special case = 1 holds for all i = 0, . . . , d according to Lemma 

3.1. Thus, the upper tail where we know the distribution of the discretized version 

V coincides with that of the initial SGPP V. 

It is obvious that the pathwise structure of the discretized version of an SMSP 
we have established in Lemma 3.3 now carries over to the SGPP case since the 
assertion in this lemma solely follows from the structure of go, ■ • ■ , and the fact 
that the initial process is nonpositive with probability one. 
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Lemma 4.4. The SGPP defined in (7') fulfills for i ^ 1, . . . ,d 



and 



sup Vt = max (l^.^i , T4, ) 

te[si_i,si] 



inf Vt = -\\{Vs^_,,Vs^)\\^ 



This minimum is attained for t ~ {si^iVs-_-^ + SiVs-)/Vs-_-i^ + Vs^)- 

Now consider a sequence of discretized versions y('') of an SGPP V with grid 
Qd, where the fineness of Qd converges to zero. Repeating the arguments in the 
proof of Theorem 4.5 yields the foUowing result. 

Theorem 4.5. The sequence of processes V'-'^\ d converges uniformly to V, 



0. 



4.2. The Mean Squared Error of the Discretized Version. The aim of this 
section is the calculation of the mean squared error of the predictor Vt of Vt . We 
obtain again some kind of pointwise convergence in mean square of a sequence of 
discretized versions with increasing fineness to the initial SGPP. Nevertheless, there 
is a difference to the considerations in the previous section. In contrast to the case 
of max-stable distributions, we typically only know the distribution of a GPD rv in 
the upper tail. Note that the function W{x) := 1 — ||a;||^, a; < 0, ||a;||^ < 1, does 
not define a multivariate df in general, see cf. Falk et al. (2010). This fact forces 
us to consider conditional expectations in this section. 

In the bivariate case, however, W defines a df, and we can assume that a GPD has 
this representation on the whole domain. The next Lemma is on some conditional 
moments of bivariate standard GPD rv in general. 

Lemma 4.6. Let ([/, V) be a standard GPD rv, i. e. there exists some D-norm 
\\-\\]j such that P(U < u,V < v) = 1 — u,v < 0, \\{u,v)\\jj < 1. Then we 

have for all such u, v 

(i) 

P{U>u,V>v) = \\{u,v)\\,~\\M\\^, 
and, in case of \\-\\^ ^ ||-||^, 
(n) 

E{U^\U >u,V >v) 

1^3 + ^2 + WMW^) + «3 ;;A' ^u/v ||(^ax(s,i), 1)11^ ds dt 
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(Hi) 

fjl\Umn ds dt + v^j;'''\\{r,l)\\^ dr 



E{UV\U > u,V > v) ^ - 



IK") 111 - III") ''J II n 

Note that the case ||-||^ = \\-\\-^ has to be treated with caution. It represents 
the case of uniform distribution on the hne {(x, y) : a;, y < 0, x + y = —1}, which 
means that no observations fall in any rectangular [it,0] x [w,0], u + v > —1, cf. 
Falk et al. (2010). 

Proof of Lemma 4-6. (i) We have 

PiU >u,V >v) = 1- P{U <u)- P{V < v) + P{U <u,V <v) 
= l-(l + u)-(l + «) + l-||(u,«)||^ 
^\\{u,v)\\,~\\{u,v)\\,,. 

(ii) We obtain by Fubini's theorem and elementary computations 

E {^{U>u,V>v]U'^) 

x^ {P*{U,V)){dix,y)) 

u,0]x[v,0] 

if I l[.fl]{s)-l[.fi]{t)dsdt\ {P*{U,V)){d{x,y)) 

[ufi]x[v.O] \Ju Ju J 

" ' I / 1m](s)-1[.,o]W (P*(C/,F))(d(x,y)) 1 dsdt 

ffil {P*{U,V)){d{x,y))\ dsdt 

'u Ju W [«,min(s,t)] X [t),0] / 
^0 



I I P(C/ G [u,min(s,t)],y e [w,0]) ds 

J u J u 

I P{U < min(s, t)) ds dt- / P {U < u) ds dt 

J u J u J u J u 

+ / P{U <u,V <v) ds dt~ / PiU <inm{s,t),V <v) ds dt 

J u J u J u J u 



rO /-O 

1- ll(5,t)|loo ds dt-u^l + u) 

' u J u 
,2 



+ u'{l-\\{u,v)\\j,)- j I l-\\{mmis,t),v)\\j, ds dt 
2 



rO rO 



3 



-u'^ {u+\\{u,v)\\^) + I I \\imm{s,t),v)\\jj ds dt 
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2 f-u/v f-u/v 

= --u^ -u^{u+\\{u,v)\\jj)-v^ I I ||(max(ri,r2),l)||^ dri dra, 
-J Jo Jo 

where wc substitute ri ~ s/v and r2 — t/v in the last equaUty. Together with 
assertion (i), the statement is proven, 
(iii) Similar arguments as in the proof of (ii) yield 

E{l{u>uy>v]UV)^ / P{U ^[u,slV e[v,t]) dsdt 

J V J u 

= 11 P(U <s,V<t)dsdt~ / P{U <s,V<v)ds dt 

- / / P{U < U,V < t) ds dt + / PiU <U,V < v) ds dt 

rO M M 

\{s,t)\\D dsdt+ / / \\{s,v)\\^ dtds 



, I.; II ^ Li., uu ^ I II ^a, L-yii jj 

V J u J u J V 

nO nO nO nO 

/ / KuMd dsdt- / 



(it, v)\\jj ds dt 



pO 



\\{s,t)\\^ dsdt-v / \\is,v)\\^ ds 



u / \\{u,t)\\D dt~-uv\\{u,v)\\jy 

J V 

' \\is,t)\\^ dsdt-v' f ^ \\ir,mn dr 

V Ju Jo 
/•■"/« 

\\{l,r)\\^ dr-uv\\{u,v)\\^. 

Jo 



□ 



Example 4.7. In case of total dependence of U and V (i.e. ||-||^ = IMIooJ and 
u = V =: c, the formulas in Lemma 4-. 6 (ii) and (ii) become 

EiU^W >u,V>.) = _M±f!(^^l±^ ^ 

— c 3 



- f f min(s, t) ds dt + c' + c' - 5 o 

E(UV U >c,V>c) = — = -c^. 

— c 3 

We now return to the discretized Version V of an SGPP V. Again, we have 
to show that for every t 6 [0, 1] the rv {Vt, Vt) follows a standard GPD. The exact 
same arguments as in Lemma 3.8 provide the bivariate df of this rv, at least in the 
upper tail. 
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Lemma 4.8. Let V — (^Ote[o,i] an SGPP with generator Z ~ {Zt)te[o,i]- 
Denote byY = (^t)te [0,1] discretized version with grid {sq, ■ ■ ■ , Sd} and generator 
Z = iZt)te[o.i]- Then {Vt,Vt) defines a bivariate standard GPD rv for every t £ 
[0, 1]. Its df is given by 

P{Vt<x,Vt<y) = l-\\{x,y)\\o, 



1 - 



[x,g,.i{t)y,g,{t)y) , t e [s,_i, s^], i ^ 1, . . . , d, 

V / Dt.^-l., 



for X, y close enough to zero, where ||-||^^ is the D-norm gernerated by {Zt, Zt) and 
ll'llj^^ . ^ . is the D-norm generated by {Zt, Zsi_-^, Zg.). 

We close this section with the calculation of the mean squared error of Vt , under 
the condition that Vt and Vt attain values that are close enough to zero, such that 
we have a representation of the df of {Vt,Vt) available in this area. Again, this 
means squared error converges to zero. 

Suppose V is an SGPP with generator Z and choose a sequence of grids of 
the interval [0, 1] with fineness converging to zero as d increases. Denote by V'''^\ 

G N, the sequence of dicretized versions of V with grid Gd, and by Z^"^', d G N, 
their generators. Denote further by ||-||^(ci) the D-norm generated by {Zt,Z^t''^), 
te [0,1], deN. 

Theorem 4.9. Let V and V'-'^\ d e N, be as above. Suppose |Ml^(d) ^ 
d G N. Then we have for c close enough to zero 



((' 



E[ [Vt-Vi^^ 



Vt > c, V}''^ > c ] ^d^oc 0. 



Proof. According to Lemma 4.8, the rv ^Vt, V'/'^''^ , c? G N, is a standard GPD rv 
with pertaining D-norm ||-||n(d). We have already shown in the proof of Theorem 
3.9 that ||-||^(d) ->d^oo IMItxj Substituting ||-|j^ by || -jlj^ in the numerators 

of Lemma 4.6 (ii) and (iii) leads to finite integrals exclusively, which is why we can 
apply the dominated convergence theorem in each of these integrals. Therefore, we 
obtain by the calculations in Example 4.7 for all t G [0, 1] 



E[(Vt- f/'*^ 



Vt > c, vy^ > c 



= E[V^\Vt > C,V}^^ > c) - 



:) - 2D (VtVt^''^\Vt>c,Vt^''^ >c) 



Vt > c, V}""^ > c 



10 o 10 o 

— c"' 



□ 
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